Smectic ordering in liquid crystal - aerosil dispersions II. Scaling analysis 
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Liquid crystals offer many unique opportunities to study various phase transitions with continuous 
symmetry in the presence of quenched random disorder (QRD). The QRD arises from the presence 
of porous solids in the form of a random gel network. Experimental and theoretical work support 
the view that for fixed (static) inclusions, quasi-long-range smectic order is destroyed for arbitrarily 
small volume fractions of the solid. However, the presence of porous solids indicates that finite-size 
effects could play some role in limiting long-range order. In an earlier work, the nematic - smectic-A 
transition region of octylcyanobiphenyl (8CB) and silica aerosils was investigated calorimetrically. 
A detailed x-ray study of this system is presented in the preceding Paper I, which indicates that 
pseudo-critical scaling behavior is observed. In the present paper, the role of finite-size scaling and 
two-scale universality aspects of the 8CB+aerosil system are presented and the dependence of the 
QRD strength on the aerosil density is discussed. 

PACS numbers: 64.70.Md,61.30.Eb,61.10.-i 



I. INTRODUCTION 

The study of the effect of quenched random disorder 
(QRD) on phase transitional behavior remains an at- 
tractive area of research due to the broad implications 
outside the laboratory. The underlying physics has ap- 
plications ranging from unique assemblies of complex flu- 
ids to doped semiconductors. Many systems have been 
the focus of both theoretical and experimental studies. 
The experimental efforts have concentrated on idealized 
model systems in the hopes of isolating the essential fea- 
tures of quenched random disorder. They include the 
still enigmatic superfluid transition of 4 He in aerogels 
and porous glasses, the superfluid transition and phase 
separation of 4 Hc- 3 Hc mixtures in silica aerog els and 
doped magnet systems Relatively recent efforts with 
liquid crystal (LC) - silica composites ||, |[ [|, [1 @' q El 
have demonstrated that these are especially interesting 
model systems. They are of particular importance as a 
way to access "soft" (elastically weak) phases of continu- 
ous symmetry, which are directly coupled to surfaces and 
external fields. 

The general consensus is that the physics of QRD in 
liquid crystals is essentially contained by a Random-Field 
approach [^0| . Recent theoretical efforts predict that 
an Ising system with quenched random fields will move 
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towards a new Random-Field Ising (RFI) fixed point 
with increasing disorder. However, a Random-Field XY 
(RFXY) system has no new fixed point that is stable. 
Here, with increasing strength of the disordering random 
field, an RFXY system still has flows toward the XY 
fixed point until long-range order is destroyed JlC| ] . Thus 
one expects, in general, that a 3D-XY system subject to 
random-field perturbations has no true long-range order 
(LRO). A detailed theoretical study of QRD effects on 
smectic ordering in liquid crystals [ [hi] concludes that ar- 
bitrarily small amounts of QRD destroy even quasi-LRO 
and hence, no true smectic phase exists since the smec- 
tic correlation length remains finite for all strengths of 
disorder and all temperatures. This theoretical conclu- 
sion is in agreement with a recent x-ray study of octyl- 
cyanobiphenyl (8CB)+aerosil dispersions |^|, the detailed 
results of which are presented in the companion paper 
to this work, denoted hereafter as paper I || . This 
x-ray study reveals a finite, though large, smectic corre- 
lation length for all temperatures and densities of silica. 
However, smectic thermal fluctuations still exist above a 
pseudo nematic to smectic transition at T* (close to but 
below T^ A for pure 8CB). These smectic fluctuations are 
expected to remain in the XY universality class but also 
show crossover behavior from Gaussian tricritical (TC) 
to 3D-XY with increasing strength of disorder fil| . 

In all fluid systems studied to date as models of such 
QRD effects, including the liquid crystal system men- 
tioned above, the random perturbations are introduced 
via the embedding of a random (gel-like) solid structure 
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into the phase ordering material. An open question re- 
mains as to the connection between the concentration of 
such solid inclusions and the strength of the random dis- 
ordering field. Also, the identification of QRD is compli- 
cated by finite-size effects which could, in principle, play 
a dominant role in such systems. In simple finite-size 
scaling (FSS), where the confining surfaces play no in- 
teractive role, the bulk critical correlation fluctuations are 
cut-off at a length dictated by the distance between sur- 
faces, which corresponds to a minimum reduced tempera- 
ture where the transition is "truncated" . However, when 
the surfaces are arranged in a random manner with high 
void connectivity in order to introduce QRD, the dis- 
tance between surfaces no longer acts as an upper length 
scale in the system, and changes in the transition's criti- 
cal behavior may also occur. Given the absence of LRO 
in such perturbed systems, the required characterization 
of the critical behavior may not be possible. In spite of 
this, if a critical power-law analysis of the transition heat 
capacity data is available, then, through two-scale uni- 
versality, the critical behavior of the correlation length 
for T > T* may be estimated and compared with direct 
measurements. Finally, if the introduced random sur- 
faces have in addition the freedom of an elastic response, 
then coupling between the gel and host elasticities can 
occur. This latter effect has only begun to be explored 
theoretically @ 0, [f§. 

The particular system that is the focus of this paper is 
a silica colloidal gel of aerosil particles dispersed in a liq- 
uid crystal denoted as LC+aerosil. The analysis will also 
be applied to another, earlier, system of an aerogel (fused 
silica gel) structure embedded within a liquid crystal de- 
noted as LC+aerogel. The two are nearly identical in 

every respect fractal-like nature of the gel structure, 

surface chemistry, and density save for their relative 

elasticity. Additionally, the ease at achieving nearly ar- 
bitrary silica densities for the aerosil system allows for 
greater control of the disorder. For the 8CB+aerosil sys- 
tem, thermal evidence for two regimes of behavior has 
been found [Q: low-density gels where pseudo-critical 
behavior is closely related to that for the pure LC and 
higher density gels where all transition features appear to 
be smeared. More rigid aerogels in LC+aerogel systems 
are crudely like the LC+aerosil gels in the high density 
regime but differ in some important ways since the elas- 
tic strain imposed by the random anchoring surfaces of 
aerogels is fully quenched. It appears that the disorder 
introduced by an aerogel is so great that all "transition" 
features are dramatically smeared and the physics of such 
systems may be more closely related to static random- 
elastic strain disorder. 

In spite of the loss of smectic LRO, Paper I [12] and 
an earlier calorimetric study of 8CB+aerosils show 
that smectic thermal fluctuations still play an impor- 
tant role in LC+aerosil systems. Many concepts from 
pure materials such as finite-size scaling, two-scale uni- 
versality, and tricritical-to-XY" crossover due to variable 
de Gennes coupling need to be considered in addition 



to intrinsic new quenched random effects that dominate 
at temperatures below an effective transition tempera- 
ture. The present paper is organized as follows. The 
relevant characteristics of an aerosil gel are described in 
Section II. Section [II reviews the essential features of 
LC+aerosil behavior near the N-SmA transition. The 
application of scaling analysis to calorimetric and x-ray 
results on the 8CB+aerosil system is presented in Sec- 
tion IV with comparisons made to the 8CB+aerogel sys- 
tem. Section |y| summarizes the conclusions that can be 
drawn from such scaling analysis and discusses the is- 
sue of the relationship between the RF strength and the 
aerosil density. Appendix |A| reviews the relevant N-SmA 
critical behavior in pure liquid crystals, and Appendix [b| 
presents the necessary background for two-scale univer- 
sality and finite-size scaling. 



II. GENERAL DESCRIPTION OF AEROSIL 
GELS 

For liquid crystal systems, the introduction of 
quenched random disorder typically requires the inclu- 
sion of random solid surfaces. This can be accomplished 
by the percolation of a low-volume-fraction gel struc- 
ture randomly arranged throughout the LC host. Such 
gels can be physically realized by a diffusion-limited- 
aggregation process which forms fractal-like structures 
having a wide distribution of void length scales. In prac- 
tical terms, the fractal-like character is limited to length 
scales much larger than the size of the basic unit and 
smaller than some macroscopic size limiting the gel, i.e., 
the sample size. 

Because of the hydroxyl groups on the surface of 
the 70 A diameter hydrophilic aerosil (Si02) spheres 
used in this work, hydrogen-bonding is possible between 
aerosil particles p6[ . When dispersed in an organic liq- 
uid medium, aerosil particles comprising 3 to 4 lightly 
fused spheres and having a mean radius of gyration of 
« 240 A [Q will attach to each other and form a gel 
by a diffusion- limited aggregation process. This gel can 
be thought of as a randomly crossing "pearl necklace" 
of silica and a cartoon depiction is given in Fig. The 
hydrogen-bonded nature of the silica " links" is relatively 
weak and gives these gels the ability to break easily and 
reform on moderate time scales (such gels are termed 
thixotropic). In addition, because of the diffusion- limited 
aggregation process by which gelation occurs, the struc- 
ture of the final gel may become anisotropic if the gela- 
tion occurs in an anisotropic fluid, i.e., gelation in a well 
aligned nematic or smectic liquid crystal. This gives such 
colloidal gels very attractive uses in future research as a 
route to studying anisotropic random disorder p7| . 

These aerosil gels are very similar in structure to the 
well-known and previously studied aerogels, which are 
another type of fractal silica gel. Aerogels are formed 
by a reaction-limited aggregation process and form a gel 
nearly identical to that of aerosils except that the ba- 
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FIG. 1: Circles and "hairs" (upper left) represent type-300, 
70 A diameter, aerosil particles and 8CB molecules, respec- 
tively, drawn to approximate scale. This cartoon corresponds 
to an average void length l a ~ 400 A and ps ~ 0.20, where 
the density units are grams of silica per cm 3 of 8CB. The solid 
volume fraction is $ « 0.08. Open and shaded parts of arrow 
depict void and solid chords respectively. 



sic silica units in aerogels are chemically fused together. 
Thus, aerogels have a large shear modulus (they break 
before yielding) while in contrast aerosil gels have a quite 
small, density-dependent, shear modulus [ p"8[ . Thus, 
aerosil gels can respond elastically to strains, which may 
turn out to be a crucial difference as will be discussed 
later. For the present work, the thixotropic character en- 
sures that the silica strands mainly dictate, at low silica 
concentrations, the local nematic director without impos- 
ing high-energy elastic strains over the material through- 
out the void. As the silica density increases, the aerosil 
gel eventually becomes stiff enough to elastically strain 
the host fluid. 

For any random gel structure, a mean distance be- 
tween solid (gel) surfaces or a mean void size, l , can 
be uniquely defined despite the wide distribution of void 
sizes. The definition of l in terms of macroscopic, and 
experimentally accessible, quantities begins by imagining 
a "straw" of uniform cross-section A sent through the gel; 
see Fig. The places where the gel randomly intersects 
the "straw" defines a solid length while the distance be- 
tween intersections defines a void length. The relevant 
macroscopic quantities are the specific surface area a in 
total surface area per mass (given as 300 m 2 g _1 for 
the type R300 aerosil used in the 8CB+aerosil samples) 
and the reduced density, ps = mass of solid per open vol- 
ume or in our case grams of silica per cm 3 of LC. Since 
the cross-section of this imaginary straw is uniform, the 
proper summing of the total void and solid volumes of 
the straw depends only on the void and solid lengths re- 
spectively, noting that each void must be bounded by 
two walls. The sum of the total solid length and the to- 
tal void length is simply the length of the straw which 
spans the sample. 

The two requirements for volume and length defined 



above allow the definition of the average void length as 

l = 2/aps . (I) 

See Ref. Q] and for a more detailed derivation. 
Strictly speaking, this definition of l is valid only in the 
dilute regime where the addition of more solid does not 
significantly change the specific surface area. However, 
as the concentration of solid increases, more surface area 
is lost due to multiple connections (" clumping" of basic 
units); thus the specific surface area should be a decreas- 
ing function of the solid volume fraction A limiting 
case of interest is when completely enclosed pores occur. 
Here, the open volume is now bounded on all sides and 
a re-derivation of the mean void (now pore) size along 
the lines described above yields l Q = 6/aps- This can be 
recast into the form given in Eq. (Q) if a pore — > a VO id/i. 
Thus, Eq. (JTJ) is quite general for a random gel structure 
if the variation of a with ps is known. The exact varia- 
tion of a with solid concentration depends on the specific 
process of densification of the gel. Since the gel struc- 
tures for both aerosil gels and aerogels are nearly iden- 
tical, results for void sizes determined from small-angle 
x-ray scattering (SAXS) of various density aerogels |2G] 
can be used to estimate the variation of a(ps), in n? 
per gram, as a = 300 — 103. 8ps, where ps has the units 
grams of SiC>2 per cm 3 of LC. The condition for closed 
pores is crudely predicted to occur at ps ~ 1.97 g cm~ 3 
(assuming a continuous process of diffusion-limited ag- 
gregation). For all LC+aerosil and LC+aerogel samples 
studied to date, the density of silica employed has been 
well below that limit and the use of this estimate of a(ps) 
in Eq. (Q) reproduces quite closely the SAXS measured 
void sizes of aerogels. This provides some confidence that 
a useful representation of a characteristic length scale l 
for fractal-like gels is available. 

From the above discussion, the reduced density ps is an 
attractive quantity to describe the gel and its disordering 
character pM. Note that the volume fraction of LC in 
a silica boundary layer of thickness lb is given by p = 
haps @- This quantity p is the fraction of LC filling the 
voids that is in direct contact with the solid surfaces and 
thus considered strongly "pinned". Since p is a natural 
measure of the "disordering strength" of the gel, ps is 
expected to be linearly related to this QRD strength for 
low to moderate ps values |p2| . 

III. NEMATIC - SMECTIC-A BEHAVIOR FOR 
LIQUID CRYSTAL - SILICA DISPERSIONS 

The essential feature of importance here is the obser- 
vation that low silica density LC+aerosil samples exhibit 
pseudo-critical behavior that is parallel to the critical be- 
havior exhibited by pure LCs in spite of the absence of 
smectic long range order in LC+aerosils. fn Paper 1, the 
usual scaling concepts are shown to hold for the relation- 
ship between the normalized thermal fluctuation ampli- 
tude <Ji and the parallel correlation length £|| for T > T*, 



where T* is an effective N-SmA transition temperature. 
These concepts also hold for the temperature dependence 
of the normalized integrated area of the static (QRD) 
fluctuation term in the x-ray structure factor for T < 
T*, where this area is proportional to , and we drop 
hereafter the superscript N denoting normalization. In 
the latter case, it is shown in Paper I that a2 ~ (T*— T) x , 
where the " critical" exponent x is essentially the same as 
2/3 for the smectic order parameter squared in pure liquid 
crystals. Calorimetric data for 8CB+aerosils j| support 
this view that effective critical behavior occurs for low 
silica density LC+aerosil samples. 

The background given in App. ^ for trends in N-SmA 
critical behavior for pure liquid crystals as a function 
of the McMillan ration Rm =Tna/Tni is pertinent 
to LC+aerosil systems. It appears that increasing the 
density ps of the thixotropic gel of aerosils decreases 
the smectic-nematic coupling described in App. |A|. In- 
deed, there is a simple empirical connection between 
the variable Rm for pure LCs and the variable ps for 
LC+aerosils. As described in Section IV of Paper I, an 
effective McMillan ratio R e // (siV) for LC+aerosil systems 
can be defined as 



R% f (si\) = 0.977 ~0.47p s 



(2) 



Figure 11 in Paper I demonstrates this equivalence of Rm 
and ps as measures of changes in the smectic-nematic 
coupling in pure LCs and LC+aerosils. This connection 
will be discussed further in Sec. [IvjB. 

The full power-law form in terms of the reduced tem- 
perature t =|T— T* I / T* used to analyze experimental 
specific heat data associated with the N-SmA phase tran- 
sition is If23| 



ACp(NA) = Cp(observed) 



A ± t- a {l 



C p (background) 
+ D ± t Al ) + B c , (3) 



where the critical behavior as a function of reduced tem- 
perature t is characterized by an exponent a, an ampli- 
tude A^ above and below the transition, a critical back- 
ground term B c , and corrections-to-scaling terms char- 
acterized by the coefficients and exponent Ai ~ 0.5. 
The excess N-SmA specific heat (heat capacity per gram 
of LC) data for 8CB+aerosil samples Q is reproduced in 
Fig. @ in log- log form in order to illustrate the quantities 
described above. Figure [2] highlights the quality of the 
fits to the standard power-law form for 8CB+aerosil sam- 
ples with low ps- Note the importance of the inclusion 
of the corrections-to-scaling term, seen as curvature at 
high t. In fact the quality of the fit for the lowest density 
8CB+aerosil sample rivals that observed for most pure 
LCs. In addition, the changing shape of the N-SmA heat 
capacity peak with ps is clearly evident. As discussed 
in Ref. E|, it is clear from Fig. || that a fit to the data 
with Eq. (|3|) can include only data for | t > since 
the AC P peak is truncated at a finite maximum value 
h M = AC™ ax {NA). 
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FIG. 2: Specific heat due to the N-SmA phase transition, 
ACp(NA), as a function of reduced temperature t, for three 
8CB+aerosil samples with aerosil densities ps = 0.022, 0.052, 
and 0.092. Data taken from Ref. Open circles and dashed 
lines represent data and fit for t > (above transition) while 
filled circles and solid lines represent data and fit for t < 
(below transition). Fits were made using Eq. (^). Also 
indicated for each sample are hu = AC'™ ax (NA) and the 
minimum reduced temperatures for which the data can be 
fit with a power-law. 



The role of pseudo-critical behavior for the 
8CB+aerosil system is fully described in Paper I 
and Ref. B. In the present paper, the applicability 
of two-scale universality and finite-size scaling in de- 
scribing the pseudo-critical behavior in LC+aerosils 
is investigated. Essential scaling background material 
is introduced in App. |b| and these concepts are im- 
plemented in Sec. IV for the analysis of the N-SmA 
"transition region" for 8CB+aerosils. 



IV. ANALYSIS OF THE NEMATIC - 
SMECTIC-A TRANSITION IN THE 8CB+SIL 
SYSTEM 

A. Specific heat behavior 

In order to utilize the scaling analysis described in 
App. [b|, the maximum length scale £m and the appropri- 
ate critical flu ctuatio n p arameters must be substituted 
into Eqs. (B2) and (B4), which are repeated here for 
convenience 
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FIG. 3: Truncation plot of h M = ACTf{NA) versus p s for 
8CB+aerosils |j and 8CB+aerogels Q. Both the dashed 
and solid lines are FSS predictions based on using bulk 8CB 
critical parameters £ii and vn to find the reduced tempera- 
ture where £m = lo- The dashed line represents the simple 
FSS result based on using only the leading singularity in 
Eq. (^); see text. The solid line denotes FSS using the full 
heat capacity function. The xs are based on a RFS analysis 
using two-scale predicted critical parameters (see text) and 
£,m = Cf T found in Paper I. 



h M = A ± (a I /^o) a/ ^(l + D ± (CM/C\io)- Al/ ^)+B c , 

(5) 

to describe the fractional temperature rounding of the 
transition and the specific heat maximum, respectively. 
In this paper, two approximations will be explored. In 
the first case, we use the mean void size as the cut- 
off length scale £m — lo and the bulk critical param- 
eters are used. This approach represents conventional 
finite-size scaling, denoted as FSS, where the cut-off 
length scale is set by a natural length of the "confine- 
ment" . In the second case, we use the calorimetrically de- 
termined 8CB+aerosil critical parameters and two-scale 
predictions for the bare correlation length and exponent 
and set the cutoff length scale to the saturated paral- 
lel smectic correlation length found by the x-ray analy- 
sis in Paper I, i.e., £m = £,^ T ■ This analysis recognizes 
that the random-field effects truncate the growth of or- 
der and tests whether two-scale universality is obeyed 
on approaching this truncation. We label this approach 
random- field scaling, or RFS for short. 

The specific heat maximum hu = AC™ ax (NA) 
is plotted versus ps in Fig. ||| for 8CB+aerosil and 
8CB+aerogel samples. Given as a dashed line on this 
log-log plot is the simple scaling prediction using only 
the leading singularity of the pure 8CB heat capacity and 
the mean- void size as the truncation length. This ignores 
the critical background and corrections-to-scaling terms 



FIG. 4: Dependence on ps of the fractional round-off region 
(gap about T*) where C p power laws fail. The solid line de- 
picts FSS predictions using bulk critical parameters and the 
fluctuation cut-off length, £a/ = lo, while the xs depict RFS 
predictions using two-scale critical parameters and £m = £|f T 
from Paper I. 



in Eq. (||) and yields a straight line having a slope of 
a/v\\ = 0.45. This is the usual theoretical FSS predic- 
tion for hu , and its failure highlights the importance of 
using the full expression given by Eq. (|^). Surprisingly, 
the full analysis denoted as FSS, which uses the pure 
8CB critical parameters and £m = lo from Eq. ([!]), ap- 
pears to work very well over the entire range of p$. This is 
surprising since the critical behavior of the 8CB+aerosil 
samples is changing with ps as seen in Fig. and the 
saturated parallel smectic correlation length £y is much 
larger than l a for all ps ]T^ |. Thus, we are suspicious that 
this agreement may be accidental. 

The analysis denoted as RFS uses the evolving specific 
heat critical behavior and two-scale universality predic- 
tions (see App. |b|) for the equivalent correlation length 
critical behavior at each ps- The low-temperature exper- 
imentally measured saturated parallel correlation length 
£k T is used as the truncation length This analysis 
reproduces very closely the observed heat capacity maxi- 
mum and is completely consistent with both the effective 
critical behavior and the maximum smectic correlation 
length. Unfortunately, this analysis is only applicable up 
to ps ~ 0.1 since a critical analysis of AC P (NA) is not 
possible for larger ps § . Note that /*m for 8CB+aerogel 
samples can not be described by either scaling methods. 
For either scaling approximation to reproduce the aero- 
gel results, a far smaller £m is required, which indicates 
that the aerogel has a much stronger disordering influ- 
ence than the aerosil at any given ps value. 

The fractional rounding of the transition ST* /T* for 
8CB+aerosil and 8CB+aerogel samples is given versus 
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p s ( g cirf 3 ) 

FIG. 5: The ps dependence of the N-SmA transition en- 
thalpy SHna = / AC p (NA)dT. As before, filled and open 
circles are data taken from [^(J and Q], respectively. The 
solid line depicts FSS using bulk 8CB critical parameters 
and £m = lo- 



ps in Fig. For samples where critical specific heat fits 
were not possible, i.e., 8CB+aerosil for p$ > 0.1 and 
all 8CB+aerogel samples, the fractional rounding is esti- 
mated ad hoc as s=s 10% larger than the width between 
inflection points in AC p . As was seen in Fig. ||, the 
FSS analysis works well for all densities of 8CB+aerosil 
samples. As before, this is surprising given the known 
changes in critical behavior and the fact that > l . 
The RFS analysis predicts a somewhat sharper transi- 
tion than is observed. The observed rounding is not likely 
influenced by the amplitude of temperature oscillations 
employed by the ac-calorimetric technique in Refs. Q and 
p0[ , which is on the order of 5 mK and would account 
for a fractional rounding of only ~ 10~ 5 . The estima- 
tion i~ = t+ explicit in Eq. (||) may be in question as it 
assumes that the unknown critical behavior of the cor- 
relation length below T* is the same as that above T*. 
A consequence of these arguments is that the agreement 
of FSS is likely accidental although intriguing. Again, 
the 8CB+aerogel fractional rounding is much larger than 
that for the 8CB+aerosil samples, which is an indication 
that a smaller £m is required and supports the view of 
a stronger disordering influence for the aerogel than the 
aerosil. 

Figure || presents the N-SmA transition enthalpy 
SHna versus ps for 8CB+aerosil and 8CB+aerogel sam- 
ples. Unlike Km and ST* jT* which are measures of 
truncation effects on AC P very close to the peak at T*, 
SHna = J AC p (NA)dT is sensitive to both truncation 
and changes in the shape of AC P over its entire temper- 
ature range. This is clearly seen in Fig. |^ and is dis- 
cussed in detail in Ref. 0. A finite-size scaling analysis 



for SHna proceeds by integrating the available AC P crit- 
ical form approximately ±3 K about T* (\t\ ps 10~ 2 ) to 
the point corresponding to i^. As an approximation, a 
linear evolution of AC P between t+ and i~ is assumed. 
The FSS analysis, given by the solid line in Fig. ||, does 
not agree well with the data for either system. The fail- 
ure of this model is not surprising since it ignores any 
changes in the AC P critical parameters such as A^ 1 and 
a with ps- A RFS analysis would not be very meaning- 
ful for SHna since for ps < 0.1, the input parameters 
for this model automatically insure perfect agreement. 
Also, the necessary input critical parameters cannot be 
obtained for 8CB+aerosil samples with ps > 0.1 or for 
any of the 8CB+aerogel samples. 



B. Smectic-A x-ray scattering for T > T* 

Two-scale universality is reviewed in App. |B], where it 
is assumed that the 3D-XY result 

aA + (^lo) - 0-647 (6) 

should hold for 8CB+aerosil samples. To proceed fur- 
ther, two additional assumptions are adopted that are in- 
herent in the x-ray analysis presented in Paper I: [y\\ —v±) 
and (£i| /£_Lo) for 8CB+sil samples have the pure 8CB 
values of 0.16 and 2.22 respectively for all ps- For those 
8CB+aerosil samples where critical AC P (NA) fits are 
available, the three above assumptions allow the pre- 
diction of £ii , £j_o) and Vx- These parameters will 
vary somewhat with ps since the a and A + values from 
the critical heat capacity fit vary with ps- The result- 
ing £\\(ps, T) "critical" behavior for T > T* permits the 
RFS calculations of h M and ST* /T* presented Sec. |y|A. 

Given two-scale universality and the assumptions out- 
lined above, both the bare smectic correlation length and 
the effective critical exponent may be estimated from the 
critical analysis of the heat capacity data. Thus, the self- 
consistent critical behavior predicted here can be directly 
compared to those measured in Paper I. 

Presented in Fig. ^| are the experimentally measured 
£|| values as a function of reduced temperature for the 
8CB+aerosil data taken from Paper I jl2| . The two lines 
given in Fig. || show the observed critical behavior in pure 
8CB pij , where T c is used for T*, and the two-scale pre- 
diction for £|| for the 8CB+aerosil sample where the pre- 
dicted correlation lengths differ most from those of 8CB 
(ps = 0.092). Note that although i/y and £n both vary 
with ps, the predicted overall trends of £ii(t) = £\\ t~ v U 
values differ only slightly from pure 8CB. The experimen- 
tal £,\\{t) data for various ps agree well with each other 
within the scatter but all 8CB+aerosil values are consis- 
tently larger than scaling prediction. The significance of 
this observation is not known. 

As an analog to Figure 11 of Paper I, the experimen- 
tally measured heat capacity exponents a for the N-SmA 
transition in pure LCs and the effective exponents O-eff 
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t =(T-T*)/T* 

FIG. 6: Parallel smectic correlation lengths £y measured 
in Paper I as a function of reduced temperature t for 
8CB+aerosil samples with densities ps given in g cm -3 shown 
in the inset. The dashed line represent pure 8CB behavior 
while the solid line represents the most extreme 2-scale 
predicted behavior for a 8CB+aerosil sample (see text). 
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FIG. 7: Effective heat capacity critical exponent a obtained 
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values are denoted by the horizontal dashed lines. Pure liquid 
crystal values, obtained from Ref. [E3[, are plotted versus the 
McMillan ratio, R M = Tna/Tni (filled circles). 



for the N-SmA pseudo-transition for 8CB+aerosil sam- 
ples can be plotted versus Rm and ps respectively, and 
this is shown in Fig. ^. In addition, the variation of the 
mean correlation length exponent v = (z/y + 2i/±)/3 for 
pure LCs (versus Rm) and those predicted by our two- 
scale analysis (versus p$) are also plotted. Going from 
left to right in Fig. ^ corresponds to decreasing nematic- 
smectic coupling along the lines given in Ref. |23fl for pure 
LCs. The two-scale predicted critical evolution of the 
average correlation length exponent with respect to ps 
is in good agreement with the corresponding evolution 
in pure LCs with respect to Rm- This match is com- 
pletely consistent with the correlation between ps and 
Rm seen from the experimentally measured exponents a 
and x w 2(3 in Paper I and given previously in Eq. (|^). 
Fig. [7] demonstrates that the two-scale predicted behav- 
ior given in this analysis for the correlation lengths is 
consistent with other observed trends in the evolution of 
pseudo-critical behavior with ps- 



V. DISCUSSION AND CONCLUSIONS 

Despite the loss of long-range smectic order, quasi- 
critical thermal fluctuations remain important at high 
temperatures for low silica density 8CB+aerosil samples. 
In addition, two-scale universality analysis provides a link 
between the SmA quasi-critical behavior of the heat ca- 
pacity and the correlation lengths. The smectic fluctu- 
ations are modified from the pure 8CB behavior due to 
the effects of quenched-random-disorder. A collection of 
effective critical exponents for the 8CB+aerosil system 
and selected pure LCs as a function of ps and Rm, re- 
spectively, is presented in Table |. As shown in Paper I 
and here, the density ps of an aerosil gel is directly cor- 
related to the McMillan ratio Rm of pure LCs, both of 
which are indicators of the strength of smectic- nematic 
coupling. The flow of the effective critical behavior for 
the N-SmA transition shown in Table | and Fig. [7] (see 
also Fig. 11 of Paper I) as a function of this QRD induced 
decoupling is consistent with theoretical predictions that 
no new fixed point is present for the RFXY model p0| . 
For the liquid crystal - silica dispersion system studied 
here, the flow is from near a Gaussian tricritical point 
toward the 3D-XY fixed point. 

This crossover behavior is explained for pure LCs by a 
decrease in nematic - smectic coupling as Rm decreases. 
For 8CB+aerosils, an increase in p s appears to have the 
same effect. Recent work has found that aerosil gels ex- 
hibit dynamics which can couple to a host liquid crys- 
tal presumably through direct coupling to director 
fluctuations. Also, recent deuterium NMR studies of 
8CB+aerosils found no appreciable change in the mag- 
nitude of orientational order above T* for ps < 0.1. The 
reason appears to be that the aerosil particles form a 
hydrogen-bonded thixotropic 3D gel network that pro- 
vides (a) random anchoring surfaces for 8CB molecules 
and (b) because of the flexible/fragile nature of the silica 
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TABLE I: Summary of effective critical exponents for 8CB+aerosils and selected pure LCs taken from Ref. ]23]]. The McMillan 
ratio for pure LCs is Rm =T Wa/T 'mi ■ Note that ps is given in units of grams of silica per cm 3 of liquid crystal, l is in A, 
and T* is in Kelvin. The last two columns refer to the fitting parameters for = B(T*— T) x given in Paper I. Going down 
the table from XY to Tricritical corresponds to increasing the smectic - nematic coupling. The v = (y\\ + 2v±)/3 entries in 
square brackets represent the scaling predictions for the 8CB+aerosil samples. For the columns giving 2 — r) and x values for 
8CB+aerosil samples, the values for and 2/3 = 2 — a — 7 are given in parentheses for pure liquid crystals. 
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gel, random elastic dampening of elastic (director) fluc- 
tuations in the liquid crystal. Both effects will reduce 
the nematic orientational susceptibility by suppressing 
director fluctuations. Thus, increasing ps in LC+aerosil 
samples is equivalent to decreasing Rm in a pure liquid 
crystal, which in the case of 8CB drives its critical behav- 
ior towards XY. This may have important consequences 
for the bulk N-SmA behavior as theoretical efforts have 
mostly concentrated on the de Gennes type of smectic 
coupling to the magnitude of nematic order parameter, 
and it appears that the coupling to director fluctuations 
may play an important role in the crossover behavior. 

Because the random disorder is introduced by the in- 
clusion of network gel structures within the liquid crys- 
tal, finite-size effects can exist and may play a role in 
truncating thermally driven fluctuations. Such effects 
would explain the increasing suppression of the heat ca- 
pacity peak with increasing ps, which corresponds to de- 
creasing the mean distance between solid surfaces. Scal- 
ing analysis provides a good description of the maxi- 
mum heat capacity and the fractional rounding (or trun- 
cation) of the transition for all 8CB+aerosil samples. 
However, FSS analysis does not provide a good pre- 
diction of the trend with ps for the transition enthalpy 
SH NA = f AC p (NA)dT. The reason for this is the fact 
that the trend in SHna is dominated not by the trunca- 
tion of the ACp(NA) peak but by the changes in shape 
and size of AC P (NA) over its entire range, and the latter 
effect is due to crossover rather than finite-size. 




p s (gcrrf ) 

FIG. 8: Saturated 8CB+sil parallel £|| (filled circles), perpen- 
dicular £x (open circles), and mean £ = (^yCx) 1 ^ 3 (open tri- 
angles) smectic correlation lengths £ iT for T <C T* versus ps; 
data taken from Paper I. The low-temperature isotropic smec- 
tic correlation lengths £ reported in Ref. |)| for 8CB+aerogel 
samples are shown by the x symbols. The solid line is the 
estimated maximum nematic director correlation length £at 
for 6CB+aerosil [^], and the dashed line is the mean void size 
lo estimate given by Eq. (|l|). 
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Below the pseudo-transition temperature T*, the cor- 
relation lengths and the amplitudes of the thermal term 
in the smectic structure factor for 8CB+aerosils sat- 
urate and are approximately temperature independent 
|Q. Quenched random disorder imposed by the aerosil 
gel network dominates the smectic fluctuations below 
the pseudo-transition. Plotted in Fig. p| are the low- 
temperature parallel correlation lengths £^ T taken from 
Paper I plus the corresponding perpendicular values 
and the mean correlation lengths £ = (^n^) 1 / 3 . Also 
plotted are the "isotropic" smectic correlation lengths 
reported at low temperature for 8CB+aerogel samples 
, the mean- void size l based on Eq. (Q), and an esti- 
mate of the saturated nematic director correlation length 
£at measured in 6CB+aerosil samples [||. As discussed 
in Paper I, the parallel correlation length is much larger 
than l Q but smaller than £at for all 8CB+aerosil sam- 
ples studied. The first fact indicates than the smectic 
domains span many "voids" and it is thus reasonable 
to expect that their influence is of a random-field type, 
while the latter fact is physically reasonable since the 
observed smectic domains cannot be larger than the size 
of a nematic domain. However, 6CB does not exhibit a 
smectic phase, so it remains unknown how the director 
correlation length would change, if at all, due to the on- 
set of smectic order. Interestingly, the low-temperature 
8CB+aerogel correlation lengths appear to agree fairly 
well with l Q , while the aerogel AC P (NA) peak is severely 
rounded, indicating that 8CB is very strongly perturbed 
by the rigid fused silica gel. 



In spite of the applicability of the pseudo-critical scal- 
ing ideas discussed here and in Paper I, quenched ran- 
dom disorder plays a dominate role in the behavior 
of 8CB+aerosils for temperatures below T*. In par- 
ticular, consider the x-ray structure factor below T*, 
which is dominated by the " Lorentzian-squared" type 
term expected for random- field disordered systems p2| . 
Random-field theories postulate a random-field at each 
i th "spin" site of strength hi whose average is (hi) = 
but whose square defines the variance of the disorder 
A = (hi ■ hi) = \h 2 \ ]27[ ]. The x-ray analysis pre- 
sented in Paper I uses a structure factor analogous to 
that used for the analysis of random-field magnets given 
in Ref. [^tJ; see Eq. (1) in Paper I, Thus, the saturated 
(low-temperature, denoted by the superscript LT) values 
of the mean smectic correlation length £ LT = (Cn^l) 1 ^ 3 , 
and the amplitude of the thermal contribution o\ T , as 
well as the amplitude of the QRD contribution a 2 to the 
structure factor can be compared to the predicted scaling 
for random-field disorder. Note that a 2 is temperature 
dependent over the entire temperature range, thus we 
use a value for a 2 (AT) at AT = T - T* = -6 K in 
what follows. The relevant scaling relations predicted for 



random-field systems at low temperatures are 

f ~ A -,,A 

cn-e 3 (7) 

a 2 ~ £° = constant 
a 2 /a l ~ <T 3 

where v& = l/(d c — d), d c is the lower critical dimension, 
and d = 3 is the physical dimension of the random-field 
system ^7j. We have dropped the superscript LT for 
convenience here and in the rest of this discussion. For a 
pure, continuous symmetry, XY system that exhibit true 
long-range order d c (XY) = 2, which shifts upward for 
a random-field XY system to d c (RFXY) = 4 |§. Al- 
though the N-SmA phase transition is not a simple mem- 
ber of the 3D-XY universality class, Eqs.(Q) are expected 
to be reasonably applicable for smectics. Thus v& = 1 for 
the divergence of the smectic correlation length with the 
strength of the random-field. Substituting this value of 
va and eliminating £ from the expressions for the scatter- 
ing amplitudes in Eqs. (0), we find the predicted scaling 
relations solely in terms of the random-field variance for 
an RFXY system as 

e-A-^H- 2 

o~i ~ A- 3 ~ \h\- e (8) 
ci2 ~ A° = constant 
a 2 /a 1 ~ A 3 ~ \h\ 6 

The scaling of these quantities with respect to ps can 
be compared to the predicted scaling behavior with re- 
spect to A in order to make a connection between A 
and ultimately \h\ with ps- In Paper I the simple as- 
sumption that A ~ ps is made, while the assumption 
ps ~ \h\ ~ A 1//2 has been made elsewhere Q. 

As discussed in Paper I in terms of £|| , a one-regime ps 
analysis (a single trend over the entire range of ps values 
studied) yields for the 8CB+aerosil system £ ~ P5 10±0 ' 2 . 
Similarly, the scattering amplitude ratio as a function of 
£ yields a 2 ja x ~ £-2-68±o.26 at AT = _ 6 R The gcatter 

in the experimental values of 0,2 (AT) as a function of ps 
only allows a very rough estimate of its scaling, but it ap- 
pears to be very weakly dependent on either correlation 
length or silica density. Unfortunately, the uncertain- 
ties in o~\ T values arising from fitting the x-rays profiles 
at low temperatures (where the quenched random term 
greatly overshadows the thermal fluctuation term) and 
the added uncertainties from normalization make it un- 
suitable for testing as a function of ps- However, the 
ratio a 2 (lxT) I o-\ T is better characterized since it is in- 
dependent of the normalization procedure. 

One possible explanation for the observed systematic 
deviations from the simple power-law forms given above 
is that there is a finite maximum aerosil density above 
which no smectic correlations can exist. This maximum 
ps corresponds to a minimum mean void size, l c , below 
which the liquid crystal has insufficient space to form 
smectic layers. Since at low densities the void size l Q is 



inversely related to pg as shown in Sec. II, we propose 
an empirical relationship 

I = A(l - l c ) (9) 

for a linear scaling of the mean correlation length with 
the corrected void size l — l c . The variation of both £ and 
fl2 (T* — 6 K)/<7i with (l Q — l c ) is given in Fig. |j| where 
A ~ 0.93 and l c ~ 80 A were found by a best linear fit 
of £ with Eq. (0) J29|. The solid line in the top panel 
shows the result of this fit. In addition, the solid line 
in the lower panel of Fig. || shows the expected behavior 
of <22 (AT) /erf 1 " given the empirical relationship of Eq. || 
with l c = 80 A and the random-field scaling predicted in 
Eq. (§). A fit with a 2 (AT)/of T - (l Q - l c )v, allowing y 
to be a free parameter and fixing l c = 80 A, yielded the 
dashed line in the lower panel and the exponent value 
y = —2.6. Since this small l c value is approximately 
twice the smectic partial bilayer thickness in 8CB, it is 
reasonable that no smectic ordering can occur when l a < 
l c . The overall behavior shown in Fig. ^ as compared 
to the random-field scalings for an XY system suggests 
that the random-field strength A = \h\ 2 ~ pg, at least for 
small pg (and thus large l Q ). This use of the empirical 
relationship in Eq. @ allows us to bring the observed 
results into better agreement with predictions for scaling 
in random-field systems. 

There is a second possible explanation for the clear 
deviations from a simple power-law dependence on pg 
shown in Paper I and here in Fig. |^. This second pos- 
sibility is the existence of a more complex relationship 
between A and pg. The correlation lengths appear 
to have a weaker pg dependence than the mean void 
size below pg — 0.1, while above this density the op- 
posite occurs |S0"| . In addition, the value of the ratio 
a2 (AT = —6 K)/crf T exhibits distinctly different power- 
law dependencies on pg above and below this silica den- 
sity; see Fig. 10 and also Fig. 12 in Paper I for further 
details. Using separate power-law characterizations for 
low density (pg < 0.1) and high density (pg > 0.1), we 
find 

f ~ Ps^^ilowpg), ~ p s 1A±01 (htghp s ) (10) 
oaM ~ p s - 2±01 (lowp s ), ~ Ps 7±0 S (high Ps ) 

These two-regime fits are shown for £ and a 2 /eri in 
Fig. [l0| Recall that the concept of a low density regime 
(pg < 0.1) and a high density regime (pg > 0.1), shown 
in Eq. (|io|), is supported by the specific heat data |Q, 
where power-law fits were possible when pg < 0.1 but 
not when pg > 0.1. 

In order to establish a connection between the disorder 
variance A (a measure of the overall disordering strength 
of the individual random fields hi, A = \h\ 2 ) and the sil- 
ica density pg for low and high pg regimes, we compare 
Eqs. (§) and @ for f and the ratio a 2 (AT)/of T , which 
are the two best characterized quantities. For pg < 0.1 
the result is A ~ p s 5 . In contrast, for pg > 0.1 the result 
is A ~ pg . This idea of two regimes may be consistent 




8/ =/ o -/ c (A) 

FIG. 9: The low-temperature mean smectic correlation 
length f = (£||£±) 1//3 (t°P panel) and scattering amplitude ra- 
tio 02/01 (bottom panel) for 8CB+aerosil versus SI = (l — lc)', 
data taken from Paper I. The lower critical length-scale be- 
low which not even short-range smectic order survives is 
l c — 80 A. The solid lines depict the empirical trend £ ~ SI 
and the a 2 (AT = -6 VL)/a± T ~ Sl~ s trend expected for 
random-field scaling. The dashed line in the lower panel in- 
dicates a free fit with 02/0-1 = C(l — 80 A) v , where the value 
y = —2.6 was obtained. 

with the picture of the aerosil gel described earlier where, 
in addition to the random silica strands ("pearl-necklace" 
of aerosil beads) providing the random-field dictating the 
local orientation of the nematic director, there perhaps 
exists an elastic coupling of these tenuous strands to the 
nematic director. This coupling would dampen the size of 
director fluctuations analogous to the effect of a wider ne- 
matic temperature range, and this gives a physical inter- 
pretation to the critical flow with random-field strength 
towards the underlying XY fixed point. The apparent 
increase in the relative effect of the clastic coupling seen 
by the stronger scaling of the quenched random disorder 
strength with silica density may be an indication that the 
aerosil gel has become significantly stiffer when pg > 0.1 
suggesting the possibility that a rigidity transition has 
occurred in the gel. 

In conclusion, a combination of finite-size effects and 
two-scale universality concepts has yielded a successful 
connection between 8CB+aerosil thermal behavior and 
x-ray correlation length behavior. The truncation of 
ACp(NA) peaks, measured by h M and <5T*/T*, and the 
observation of high-temperature 8CB+aerosil correlation 
lengths that are close to those for pure 8CB are both 
well explained. The observed dependence of the inte- 
grated area SH^a on pg is not properly described by 



11 



10 3 



ltf 



H 

<5 



- 


c 




4-H 1 1 1 i i i i 

1 











0.01 0.1 1 

P s (g cm " 3 ) 

FIG. 10: The low-temperature mean smectic correlation 
length £ = (£||£l) 1//3 (top panel) and scattering amplitude 
ratio a2(AT = —6 K)/<rf T for 8CB+aerosil versus ps; data 
taken from Paper I. The solid and dashed lines depict the low 
and high density regimes, respectively, with slopes given in 
Eq. (|l0j). 



finite-size scaling, and the reason for this failure is clear. 
Trends in 5Hna{ps) reflect the changing shape of AC P 
peaks over a wide temperature range as ps varies. As 
described in Paper I [fL2| and in Sec. IV of the present pa- 



per, quenched random disorder plays a dominant role in 
changing the effective pseudo-critical exponents that de- 
scribe the smectic behavior in 8CB+aerosils. The density 
ps for 8CB+aerosil samples can be equated to the McMil- 
lan ratio Rm for pure LCs; see Eq. (||). Since in pure 
LCs a larger McMillan ratio implies a larger nematic sus- 
ceptibility, increasing Rm corresponds to increasing the 
smectic-nematic coupling. Increasing ps for LC+aerosils 
has the opposite effect, as expected. The linearity of the 
relationship between ps and Rm provides support for the 
view that the aerosil gel is involved in decoupling the ne- 
matic and smectic order parameters. It seems possible 
that the elasticity of the aerosil gel plays an important 
part in the theory of gels as a random perturbation act- 
ing on the N-SmA transition in liquid crystals. Theory 
incorporating such elastic aspects is in progress |y| . Fi- 
nally, reasonable scaling was observed of the smectic cor- 
relation length and scattering amplitudes with respect to 
the silica density or mean void size, which were roughly 
consistent with predictions for random-field type disor- 
der. However, the theoretical predictions for random- 
field systems were for uncorrelated random fields, while 
the disorder due to the aerosil gel are correlated, over 
some small length scale, due to their fractal structure 



pl| . The effect of algebraic correlations in the disorder 
on the random-field scaling predictions would be an at- 
tractive avenue for theoretical investigation. 
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APPENDIX A: RELEVANT NEMATIC - 
SMECTIC-A BEHAVIOR IN PURE LIQUID 
CRYSTALS 

1. General Description of the N-SmA phase 
transition 

Although there is still some theoretical debate, the 
best experimental evidence to date points to 3D-XY as 
the underlying critical behavior for the N-SmA transition 
over experimentally accessible ranges of reduced temper- 
ature t p^ . For pure liquid crystals, smectic ordering is 
strongly influenced by two types of coupling to nematic 
order: coupling to the magnitude of nematic order (S) 
and to the nematic director fluctuations (Sh). 



2. Smectic coupling to the nematic order 
parameter 

The first type of nematic - smectic coupling is the so- 
called de Gennes coupling of the form ip 2 S. Such a term 
affects the coefficient b of the quartic term bip 4 in the 
mean-field Landau expansion of the smectic free-energy 
in terms of the complex smectic order-parameter %p p3| . 
Since the nematic elastic constants are proportional to 
the square of S, this coupling reflects the effect of the 
"softness" of the nematic order prior to the onset of 
smectic order. The strength of this coupling depends 
on the magnitude of the nematic susceptibility, xn, and 
such a coupling can drive the N-SmA transition from 
XY-like {b > 0) to TC (b = 0) to weakly first-order 
(b < 0) with increasing xn- This is clearly seen as a 
trend with the width of the nematic temperature range, 
or the McMillan ratio Rm = Tna/^ni, which directly 
influences xn p3[ . As an example, the heat capacity crit- 
ical exponent varies from a = 0.50 to 0.10 as this ratio 
varies from Rm — 0.994 for a 2-K wide nematic range to 
~ 0.954 for a 15-K wide nematic range. For LC samples 
with large nematic ranges, Rm ~ 0.898 (45-K wide) to 
Rm ~ 0.660 (189-K wide), the experimental exponent is 
a w a XY = -0.013 pi M. 
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3. Smectic coupling to director fluctuations 

Smectic order coupling to the director fluctuations has 
the form ip 2 Sn. Thus, in addition to the "softness" of 
the nematic, fluctuations in the director orientation com- 
pete with the establishment of smectic order by excit- 
ing anisotropic elastic deformations in the smectic. The 
theory for such coupling is not yet complete nor are all 
the implications of this effect understood, but a self- 
consistent one- loop model has been put forward [ B5| . 
This type of coupling leads to anisotropy in the corre- 
lation lengths parallel and perpendicular to the nematic 
director (a feature not present in a normal 3D-XY sys- 
tem) and a very gradual crossover from a broad weakly- 
anisotropic critical correlation regime v\\ > v± (weak 
coupling limit) toward a strongly-anisotropic uu = 2v± 
regime (strong coupling limit). The strength of this cou- 
pling depends on the magnitude of the splay elastic con- 
stant Kn, which should vary as Kn ~ S 2 . Note that the 
XY model has no such splay component. Thus, a liquid 
crystal with a small nematic range will have a small Kn 
at Tna and should lie deep in the anisotropic crossover 
regime. Liquid crystals with a large nematic range will 
have a relatively large Kn at their Tna an d should 
straddle isotropic and weak-anisotropic regimes. The 
latter is observed experimentally, but strong anisotropy, 
z/|| = 2z^l, is not seen for any smectic since the narrow 
nematic range condition also induces the de Gennes cou- 
pling and thus crossover to tricritical and even first-order 
behavior |23 . 



APPENDIX B: SCALING BACKGROUND 
1. Finite-size scaling 

The concept of finite-size effects is a straightforward 
and basic idea in the modern theory of phase transitions 
p6| , [ItJ . Ignoring specific surface interactions, finite-size 
effects stem from the saturation upon cooling of the grow- 
ing correlation length in a disordered phase to a finite 
length scale. In traditional finite-size scaling (FSS), the 
maximum length is dictated by the " container" size. This 
effect truncates the transition prematurely and leads to 
three observable effects on the calorimetric data: a sup- 
pressed heat capacity maximum, a rounding of the tran- 
sition in temperature, and a suppression of the transition 
enthalpy. In addition, for a transition that breaks a con- 
tinuous symmetry, random-fields lead to a saturation in 
the growth of order. The hypothesis dictating our ap- 
proach is that the same analysis can be applied in the 
case of either type of truncation, provided changes in the 
critical behavior due to the disorder are accounted. In 
the case of random-fields, we call this analysis random- 
field scaling (RFS). 

In order to proceed, the power-law behavior of the 
smectic correlation length in the nematic phase needs 
to be considered. It is common practice to ignore 



corrections-to-scaling terms and use simple pure power 
laws although this is inconsistent with theory as dis- 
cussed in Rcf . ]32| . The correlation lengths of pure LCs 
are anisotropic with respect to the smectic layer normal 
(i.e., the nematic director for smectic- A phases) and are 
represented by effective critical exponents that are free 
parameters pi, p3] 



£ll = Z\\ot- 



(Bl) 



(B2) 



where £m and £j_ are the bare correlation lengths and 
and v±_ are the exponents parallel and perpendicular 
to the layer normal, respectively. 

For smectic liquid crystals, the parallel correlation 
length is always larger than the perpendicular, and so 
our analysis uses this length scale for the definition of the 
minimum reduce temperature. Defining the m axim um 
possible correlation length as £m> one solves Eq. (Bl) for 
the minimum reduced temperature above T* as tf n — 
(£m/£||o) -1 ^ 1/|1 ■ It is not possible to define a similar min- 
imum reduced temperature below the transition since 
the critical correlation length behavior below T* is not 
known. Thus, the equation for the fractional rounding 
(truncation) of the transition due to finite length effects 
is estimated to be 



(B3) 



6T*/T* = + |f-|) a 2t+ = 2(£ M /£|| ) 

Substituting i+ into Eq. (||) gives the relationship for the 
heat capacity maximum hu at the N-SmA transition as 

h M = A ± (£m/£|| ) qA ' 11 (1 + D^^m/^o)-^ ) + B c . 

(B4) 

Because of the importance of corrections-to-scaling for 
the analysis of AC P (NA), a log- log plot of Hm — B c ver- 
sus £m would not yield a straight line of slope ol/vh. The 
FSS effect on the N-SmA transition enthalpy is obvious 
since it involves replacing the singular AC P (NA) peak 
between i+ and t~ by hu and thus decreasing the inte- 
gral of AC P (NA) over T. 



2. Two-scale universality for T > T* 

Two-scale theory of critical phenomena relates the 
non-universal coefficient of the heat capacity's leading 
singularity to the non-universal bare correlation volume 
p5| , p8[ , and it also yields the hyperscaling relation be- 
tween the critical exponents a and v. For liquid crystals, 
there is an anisotropic version of hyperscaling fl39| 



2v\ 



(B5) 



which is empirically supported by the somewhat scat- 
tered available data on pure LCs like 8CB p3|. If such 
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hyperscaling holds, then a two-scale relation can be writ- 
ten as 

aA+(^ o e ±0 ) = k B (R+f = 13.8(i?+) 3 , (B6) 

where ks is the Boltzman constant and the value 
13.8 pertains when the correlation lengths are in units 
of A and A + is in units of J K _1 g^ 1 . The quantity 
R~£ has a different universal value for each universality 
class. For the 3D-XY model, the value of 13.8(i?^) 3 
is 0.647 [|U and several pure LCs have been shown to 



have aA + (£||o6x ) values close to this |32| : the value for 
8CB is 0.651 jl|, [24). Assuming that hyperscaling and 
XF-like pure LC values of hold for 8CB+sil samples 
independent of ps, then 

aA+(Z llo e ±0 )^ 0.647. (B7) 

Thus, heat capacity critical behavior (a and A + values) 
for LC+aerosils will allow the determination of v\\ + 2v± 

and (CuoCio)- 
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